ELSEVIER 


Available  online  at  www.sciencedirect.com 


SCIENCE 


DIRECT9 


Journal  of  Power  Sources  143  (2005)  48-56 


JOURNAL  OF 


www.elsevier.com /locate /jpowsour 


Transient  analysis  of  reactant  gas  transport  and  performance  of 

PEM  fuel  cells 

Wei-Mon  Yana’*,  Chyi-Yeou  Soongb,  Falin  Chenc,  Hsin-Sen  Chud 

a  Department  of  Mechatronic  Engineering,  Huafan  University,  Shih-Ting,  Taipei  223,  Taiwan,  ROC 
b  Department  of  Aerospace  and  System  Engineering,  Feng  Chia  University,  Seatwen,  Taichung  407,  Taiwan,  ROC 
c  Institute  of  Applied  Mechanics,  National  Taiwan  University,  Taipei  106,  Taiwan,  ROC 
d  Department  of  Mechanical  Engineering,  National  Chiao  Tung  University,  Hsin-Chu  300,  Taiwan,  ROC 

Received  28  May  2004;  received  in  revised  form  4  November  2004;  accepted  9  November  2004 

Available  online  2  February  2005 


Abstract 

The  dynamic  performance  of  PEM  fuel  cells  is  one  of  the  most  important  criteria  in  the  design  of  fuel  cells  with  application  to  mobile 
systems.  To  study  this  issue,  we  extend  our  previous  steady  model  of  gas  reactant  transport  to  an  unsteady  one  and  employ  it  to  examine 
the  transient  transport  characteristics  and  the  system  performance  of  the  PEM  fuel  cells.  With  the  assumption  of  the  two-dimensional  mass 
transport  in  the  cathode  side  of  PEMFC,  the  effects  of  the  channel  width  fraction,  X  =  £c/£b,  the  porosity  of  the  gas  diffuser  layer,  e\,  and  the 
surface  overpotential  of  the  catalyst  layer,  rj,  on  the  transient  characteristics  of  the  resultant  current  density  and  mass  transport  are  focused  in 
this  work.  It  is  disclosed  that  an  increase  in  7,  or  rj  may  lead  to  a  faster  dynamic  response  for  the  fuel  cell  when  the  PEM  fuel  cell  system  is 
started  up.  Results  of  a  typical  case  demonstrate  that,  although  the  dynamic  response  time  may  be  as  long  as  10  s  due  to  the  mass  transport 
lag,  the  fuel  cell  system  needs  only  less  than  0.4  s  to  reach  the  90%  response. 

©  2005  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

Proton  exchange  membrane  (PEM)  fuel  cells  have  been 
recognized  as  a  promising  candidate  of  power  generators 
for  both  stationary  and  automotive  applications  due  to  their 
merits  of  low  operating  temperature,  high  efficiency  and  ex¬ 
tremely  low  pollutant  emission.  In  applications  of  PEM  fuel 
cells  to  mobile  systems,  the  dynamic  response  to  the  start-up 
of  the  system  is  important.  In  addition,  the  time  interval  from 
the  start  of  operation  to  the  steady  state  is  one  of  the  signifi¬ 
cant  performance  characteristics  of  the  system.  To  attack  the 
transient  behaviors  of  the  PEM  fuel  cell,  the  time  interval  and 
other  potential  time  constants,  the  unsteady  reactant  gas  trans¬ 
port  in  the  PEM  fuel  cell  shall  be  considered.  The  dynamic 
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response  of  the  fuel  cell  depends  strongly  on  the  operating 
conditions,  transport  phenomena  in  the  cells,  electrochemi¬ 
cal  reaction  kinetics,  mechanical  design  and  manufacturing 
process,  etc.  Among  these,  the  reactant  gas  transport  in  the 
flow  distributor  and  the  gas  diffusion  layer  (GDL)  before  oc¬ 
currence  of  the  chemical  reaction  is  one  of  the  crucial  points. 

The  GDL  morphology,  geometry  and  porosity,  as  well  as 
the  flow  channel  arrangement,  e.g.,  parallel,  serpentine  or  in- 
terdigited  channel,  have  significant  influences  on  optimiza¬ 
tion  of  the  cell  performance.  In  the  past  years,  studies  related 
to  the  above  issues  have  been  performed.  Bernardi  and  Ver- 
brugge  [1]  developed  a  one-dimensional  model  to  examine 
the  proton  transport  behaviors  in  the  membrane.  Their  re¬ 
sults  indicated  that  if  the  membrane  is  maintained  at  fully 
saturated  condition,  the  resistance  then  becomes  significant 
at  the  current  densities  greater  than  200  mA  cm-2 .  Okada  and 
colleagues  [2-7]  studied  water  transport  with  and  without  im- 
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Nomenclature 

D  Mass  diffusivity,  m2  s-1 

F  Faraday  constant,  96487  C  mol-1 

I  current  density,  A  m-2 

/o  maximum  current  density  at  X  =  0,  A  m  2 

j  transfer  current  density,  A  m-2 

i b  base  width  of  the  flow  distributor,  m 
i c  gas  channel  width,  m 

i s  shoulder  width,  m 

Mo2  molecular  mass  of  oxygen,  32.0  kg  kmol-1 
N  number  of  flow  channels 

n  number  of  electrons  in  transfer  current  density 

relation 

P  pressure,  Pa 

R  universal  gas  constant,  8 . 3 1 4  J  mol- 1  K- 1 

T  temperature,  K 

t  time,  s 

5c  thickness  of  catalyst-layer,  m 

5’d  thickness  of  diffusion  layer,  m 

sm  thickness  of  membrane  layer,  m 
V  cell  voltage 

w  oxygen  mass  fraction 

wQ  oxygen  mass  fraction  at  inlet 

Greek  letters 

a  transfer  coefficient  for  the  reaction 

£  \  porosity  of  diffusion  layer 

£2  porosity  of  catalyst  layer 

£i  volume  fraction  of  ionomer  in  the  catalyst  layer 

0  surface  overpotential,  V 

k  reaction  rate  constant,  mol  m-3  s-1 

A  channel  width  fraction,  lc/l b 

Aw  membrane  water  content,  kmol 
H20/kmolS03~ 

p  density  of  oxygen,  kg  m-3 

a  ionic  conductivity  of  the  ionomer,  m-1 

x\  tortuosity  of  diffusion  layer 

T2  tortuosity  of  catalyst  layer 

Ti  tortuosity  of  the  ionomer  in  catalyst  layer 

<fi  ionomer  phase  potential,  V 


purity  ions  in  membranes  of  PEM  fuel  cells.  Gurau  et  al.  [8] 
proposed  a  one-dimensional  model  to  examine  the  character¬ 
istics  of  the  mass  transport  of  reactant  gas  in  a  half-cell,  in 
which  the  effects  of  the  porosity  and  the  tortuosity  of  GDL 
and  the  catalyst  layer  were  investigated. 

For  more  details  of  the  steady-state  flow  distribution 
and  gas  diffusion,  a  number  of  studies  on  two-  and  three- 
dimensional  transport  phenomena  in  PEM  fuel  cells  have 
been  carried  out.  Dutta  et  al.  [9]  used  the  software  package  to 
analyze  the  distributions  of  gas  density  and  velocity  in  flow 
channel  and  GDL.  Nguyen  and  coworkers  [10,11]  developed 


a  two-dimension  model  to  examine  the  water  and  thermal 
management  of  PEM  fuel  cells  of  interdigitated  flow  chan¬ 
nel.  They  found  that  the  interdigitated  gas  distributors  can 
reduce  mass-transport  overpotential  and  improve  the  flood¬ 
ing  phenomena  in  the  cathode  side.  In  addition,  it  is  found 
that  the  cell  performance  can  be  enhanced  by  decreasing 
the  electrode  thickness  or  increasing  the  number  of  the  gas 
channels.  To  obtain  an  optimization  of  the  flow  field  and 
to  improve  the  fuel  cell  performance,  Kee  et  al.  [12]  em¬ 
ployed  the  Navier-Stokes  equations  as  the  flow  model  and 
the  Michaelis-Menten  type  two-step  kinetics  as  the  gas  re¬ 
action  model  at  the  anode  to  explore  the  effects  of  flow  dis¬ 
tributors.  For  flow  distribution  in  channel  networks  of  typ¬ 
ical  planar  fuel  cell  layers  and  stacks,  Hontabon  et  al.  [13] 
proposed  a  generalized  steady- state  computational  model,  in 
which  only  mass  and  momentum  transfer  are  considered.  Um 
et  al.  [14]  studied  two-dimensional  transport  phenomenon  of 
the  PEM  fuel  cells  numerically.  The  results  showed  that  the 
hydrogen  dilution  deteriorates  the  cell  performance.  Most  re¬ 
cently,  Yan  et  al.  [15]  developed  a  two-dimensional  model  to 
investigate  gas  reactant  transport  at  various  conditions  of  flow 
channel  width  ratio  and  GDL  porosity.  They  disclosed  that 
an  increase  in  either  A  or  £  may  lead  to  a  better  cell  perfor¬ 
mance.  At  relatively  low  overpotential,  better  uniformity  in 
current  density  distribution  along  the  width  of  the  cell  can  be 
attained. 

The  previous  studies  mentioned  above  are  all  concerned 
with  the  steady- state  transport  phenomena  and  the  system 
performance  of  PEM  fuel  cells.  In  application  of  the  fuel 
cells  as  power  systems  of  automobiles,  start-up  time  is  prob¬ 
ably  a  most  significant  performance  index.  However,  in  the 
past  years,  only  a  few  analyses  focused  their  attentions  on 
the  transient  characteristics  of  PEM  fuel  cells.  Francesco  and 
Arato  [16]  investigated  the  transient  response  of  an  automo¬ 
tive  power  system  by  employing  a  model  of  fuel  cell  stack  and 
a  compressor.  Their  model  and  results  showed  that  a  proper 
range  of  the  values  of  the  air  relative  humidity  and  temper¬ 
ature  must  be  maintained  for  system  start-up  under  stable 
condition.  Yerramalla  et  al.  [17]  developed  a  mathematical 
model  for  investigating  the  dynamic  performance  of  a  PEM 
fuel  cell  with  and  without  assumptions  for  linearization.  Sig¬ 
nificant  differences  between  the  linear  and  non-linear  results 
were  found.  Without  special  emphasis  on  the  transient  behav¬ 
iors  of  the  fuel  cell  systems,  Um  and  Wang  [18]  developed  a 
three-dimensional  unsteady  computational  fuel  cell  dynam¬ 
ics  (CFCD)  model  for  analysis  of  interactions  between  mass 
transport  and  electrochemical  kinetics.  In  a  recent  investiga¬ 
tion  by  our  group,  Chen  et  al.  [19],  we  investigated  the  tran¬ 
sient  behavior  of  the  water  transport  across  the  membrane 
of  the  PEM  fuel  cell  to  seek  for  effective  control  schemes 
so  that  the  best  dynamic  performance  of  the  fuel  cell  can  be 
obtained.  It  is  found  that  both  a  larger  starting  operational 
current  density  and  a  smaller  operational  current  density  can 
lead  to  a  shorter  dynamic  response  time,  while  control  of  the 
humidification  of  the  fed  fuel  is  the  most  powerful  as  well  as 
the  most  feasible  scheme. 
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The  objective  of  this  work  is  to  study  the  effects  of  the  flow 
distributor  design  and  the  GDL  morphology  on  the  transient 
behaviors  of  reactant  gas  transport  and  dynamic  cell  perfor¬ 
mance  of  a  fuel  cell  system  subjected  to  a  start-up  operation. 
The  dynamic  response  time  for  the  system  state  changing 
from  the  start-up  operation  to  the  steady  state  will  be  exam¬ 
ined.  The  emphasis  of  the  analysis  is  placed  on  the  influences 
of  the  channel  width  characterized  by  a  fraction  A  =  lc/^b 
(=channel  width/base  plate  width),  and  the  porosity  of  the 
GDL.  By  applying  the  time- dependent  mass  transport  model 
developed  in  the  present  analysis,  the  effects  of  these  two  pa¬ 
rameters  and  the  overpotential  on  the  transient  characteristics 
of  the  mass  transport  and  the  dynamic  response  of  the  cell 
performance  are  explored. 


2.  Theory 

2.7.  Model  development 

A  schematic  diagram  of  the  cross-sectional  view  of  a  half- 
PEM  fuel  cell  model  is  shown  in  Fig.  1,  where  Sm ,  Sc,  and 
Sd  denote  thicknesses  of  membrane,  catalyst  layer,  and  GDL, 
respectively.  The  parameters  £c,  £s,  and  £b  =  l c  +  i s  stand 
for  the  widths  of  the  channel,  the  shoulder  and  the  base  plate, 
respectively.  The  area  ratio  between  flow  channel  and  shoul¬ 
der  is  one  of  the  most  important  geometry  parameters  on 
the  behaviors  of  the  PEM  fuel  cells.  With  the  assumption  of 
the  two-dimensional  mass  transport  in  the  cathode  side  of 
PEMFC,  the  unsteady  two-dimensional  transport  equations 
for  the  oxygen  mass  fraction,  w ,  in  diffuser  layer  and  catalyst 
layer  can  be  expressed  as  follows: 
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B5 
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B4 

Gas  Diffusion  Layer 
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In  the  diffuser  layer 


dw 


d 


ps  i —  =  — 
dt  dx 


sT^  pD 


dw 

dx 


In  the  catalyst  layer 


dw 


d 


ps  2 -  =  — 

dt  dx 


sTff  pD 


dw 

dx 


j 

2nF 


(1) 


(2) 


where  D  stands  for  the  mass  diffusivity  of  the  gas,  s  and 
r  are,  respectively,  the  porosity  and  tortuosity  factor  of  the 
catalyst  layer.  While  n  denotes  the  number  of  electrons  in 
transfer  current  density  relation  and  F  is  the  Faraday  constant. 
In  Eq.  (2),  j  represents  the  cathode  transfer  current  density 
and  can  be  formulated  by  the  following  equation  under  the 
assumption  that  the  oxygen  reduction  reaction  is  irreversible, 
first-order  in  oxygen  and  that  electroneutrality  holds  in  any 
representative  elementary  volume, 

j  —  uFkwqx  p 


Substituting  Eq.  (3)  into  (2)  yields 


dw 

p£25i  = 


d 

dx 


si2  pD 


dw 


+ 


d 


si2  pD 


kM0: 


dx  J  dy 
anF 


dw 

dy 


exp 


RT 


Tj  I  W 


In  the  membrane  and  the  catalyst  layer,  the  phase  potential 
equations  can  be  expressed  as  follows: 


(A)  In  the  membrane 


d 

dx 


(B)  In  the  catalyst  layer 


d 

dx 


Due  to  the  presence  of  transfer  current  density  in  the  catalyst 
layer,  there  exists  a  source  term  in  Eq.  (6).  In  addition,  S[  and 
Ti  stands  for  the  volume  fraction  and  the  tortuosity  factor  of 
the  ionomer  in  the  catalyst  layer,  and  o  for  the  membrane 
ionic  conductivity  and  is  evaluated  by 


a  =  (0.5139AW 


0.326)  exp 


1268 


(7) 


2.2.  Initial  and  boundary  conditions 


Fig.  1 .  Half-cell  model  of  a  PEM  fuel  cell.  FC  =  flow  channel;  CS  =  channel 
shoulder;  CC  =  current  collector;  BP  =  bipolar  plate;  GDL  =  gas  diffusion 
layer;  CL  =  catalyst  layer;  M  =  membrane;  B 1  ~  B7  =  boundary  numbers 
[15]. 


In  this  work,  the  fuel  cell  system  is  assumed  to  be  sta¬ 
tionary  initially  (t<  0).  At  time  t-  0,  the  fuel  cell  system  is 
subjected  to  a  start-up  operation,  thereafter,  the  system  shows 
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a  transient  response  until  the  steady  state  is  reached.  There¬ 
fore,  the  initial  condition  about  the  mass  fraction  w  is  zero.  As 
for  the  boundary  conditions  with  the  boundaries  numbered 
B1  to  B7  in  Fig.  1,  they  can  be  described  as  follows.  The 
boundary  conditions  for  the  oxygen  mass  fraction  w  are, 


w  =  wo  at  boundary  B 1 ; 
dw 

=  0  at  boundaries  B2  and  B5: 


dy 

dw 

dx 


=  0  at  boundaries  B3  and  B7: 


dw  dw 

W-  =  W- 1_,  e\lD —  =  si2  D —  at  boundary  B4. 

dy  dy 

The  boundary  conditions  for  phase  potential  0  are, 
dcp 


(8a) 

(8b) 

(8c) 

(8d) 


dy 

dcp 

dx 


=  0  at  boundary  B4: 


=  0  at  boundaries  B3  and  B7; 


0-  = 


dy 


90 

dy 


at  boundary  B5; 


+ 


0  =  0  at  boundary  B6. 


(9a) 


(9b) 


(9c) 


(9d) 


2.3.  Evaluation  of  cell  potential 


The  cell  potential  V  is  evaluated  by  subtracting  surface 
overpotential  rj  and  the  average  phase  potential  difference 
A 0  from  the  ideal  cell  potential  Voc,  viz., 

V  =  Voc  ~  0  -  A0  (10) 


and  the  value  of  A 0  can  be  evaluated  by  the  following  ex¬ 
pression, 


/o”(0 6  -  04)  d* 

4 


(11) 


in  which  04  and  06,  respectively,  denote  the  phase  potentials 
at  the  boundaries  B4  and  B6. 


2.4.  Evaluation  of  time  constant 

In  the  study  of  the  transient  behaviors  of  the  PEM  fuel 
cells,  the  time  interval  and  the  time  constant  are  important 
when  the  system  is  subjected  to  a  start-up  operation.  In  this 
work,  the  time  constant  is  defined  as  the  time  needed  for  the 
value  of  current  density  reaching  90%  of  that  at  the  steady 
state. 


3.  Numerical  method 

The  above  set  of  equations  for  time-dependent  transport  in 
PEM  fuel  cells  have  to  be  solved  numerically.  In  the  present 
study,  we  adopt  a  fully  implicit  finite  difference  method  in 


which  the  time  derivative  term  is  approximated  by  the  back¬ 
ward  difference  method  and  the  space  derivative  terms  by  the 
central  difference.  With  the  coupled  finite-difference  equa¬ 
tions,  the  linear  algebraic  equations  are  solved  by  the  Thomas 
algorithm  [20]  with  iterations.  In  this  work,  uniformly  dis¬ 
tributed  grids  having  201  and  221  points  are  used  in  the  x 
and  y  directions,  respectively.  To  obtain  enhanced  accuracy 
in  the  numerical  computations,  grid  independence  is  tested 
in  the  separate  numerical  runs.  It  is  found  that  the  deviations 
between  grid  points  of  201  x  201  and  401  x  401  with  first 
time  interval  0.001  are  less  than  2%.  Accordingly,  the  com¬ 
putations  on  the  grid  system  of  201  x  201  points  in  the  x  and 
y  directions  with  first  time  interval  0.001  are  sufficient  to  un¬ 
derstand  the  transient  behaviors  of  the  reactant  gas  transport 
in  the  PEM  fuel  cells.  To  further  check  the  adequacy  of  the  nu¬ 
merical  scheme,  the  limiting  case  of  steady  one-dimensional 
model  is  first  evaluated.  The  predicted  results  agree  well  with 
those  of  Gurau  et  al.  [8].  Through  the  above  tests,  it  confirms 
that  the  presented  model  and  numerical  method  both  are  ap¬ 
propriate  for  the  present  problem. 


4.  Results  and  discussion 

The  values  of  the  parameters  considered  in  the  present 
study  are  listed  in  Table  1.  The  numerical  results  are  pre¬ 
sented  to  discuss  the  effects  of  various  parameters.  The  ge¬ 
ometry  is  characterized  by  a  parameter  named  channel  width 
fraction, A  =  i c/t b,  and  three  values,  A  =  0.25,  0.5,  and  0.75, 
are  considered.  In  addition,  the  effects  of  the  GDL  morphol- 


Table  1 


Parameters  used  in  the  present  analysis 


Parameters 

Notation 

Value 

Temperature,  K 

T 

333 

Pressure,  Pa 

P 

303975 

Reaction  rate  constant,  molm-3  s_1 

K 

1.5  x  lO"5 

Oxygen  mass  fraction  at  inlet 

wo 

0.17 

Gas  channel  width,  m 

4 

5  x  1CT4 

Channel  shoulder  width,  m 

4 

5  x  icr4 

Thickness  of  diffusion  layer,  m 

Sd 

5  x  1CT4 

Thickness  of  catalyst  layer,  m 

V 

5  x  icr5 

Thickness  of  membrane,  m 

Vn 

2.0  X  10“4 

Porosity  of  GDL 

et 

0.8 

Tortuosity  of  GDL 

U 

1.5 

Porosity  of  catalyst  layer 

£2 

0.3 

Tortuosity  of  catalyst  layer 

T2 

1.5 

Volume  fraction  of  ionomer  in  CL 

£l 

0.3 

Ionomer  tortuosity  in  CL 

1.5 

Faraday  constant,  cmol-1 

F 

96487 

Universal  gas  constant,  kJ  kmol~ 1  K~ 1 

R 

8.314 

or  J  mol  1  K  1 

Molecular  mass  of  oxygen,  kg  kmol  1 

M 

32.0 

Number  of  electrons  in  relation 

n 

2 

participating 

Current  density,  A  m-2 

I 

Surface  overpotential,  V 

Ionomer  phase  potential,  V 

0 
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Fig.  2.  Oxygen  mass  fraction  contours  in  diffusion  and  catalyst  layers  with 
A  =  0.50  and  77  =  0.2  at  different  transients:  (a)  f  =  0.01;  (b)  7  =  0.05;  (c) 
7  =  8.35. 
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Fig.  3.  Oxygen  mass  fraction  contours  in  diffusion  and  catalyst  layers  with 
A  =  0.5  and  ??  =  0.225  at  different  transients:  (a)  7  =  0.01;  (b)  7  =  0.05;  (c) 
7  =  4.215. 


ogy  are  explored  by  considering  the  values  of  the  GDL  poros¬ 
ity  in  the  range  of  0.3  <  s  <  0.6. 

Fig.  2  presents  the  contours  of  the  oxygen  mass  frac¬ 
tion  in  the  GDL  (0<y<5xl0-4)  and  catalyst  layer 
(5  x  10-4  <  y  <  5.5  x  10-4)  with  s\  =  0.5,  77  -  0.2  and  A  =  0.5 
at  time  instants  t-  0.01,  0.05  and  8.35.  The  right  half  of  the 
x-axis  in  each  contour  corresponds  the  axial  location  of  the 
shoulder  (or  rib)  of  the  channel.  An  overall  inspection  of 
Fig.  2  discloses  that  in  the  region  corresponding  to  the  chan¬ 
nel  central  region  (left  portion  of  the  contour),  the  relatively 
higher  oxygen  mass  fraction  appears  and,  then,  the  oxygen 
diffuses  towards  the  porous  material  covering  the  rib  where 
there  is  no  oxygen  source  on  the  interface  y  =  0.  In  the  early 
initial  stage,  see  Fig.  2(a),  only  a  higher  oxygen  mass  frac¬ 
tion  presents  in  the  region  near  the  inlet  of  channel  central 
region  only.  This  is  because  the  oxygen  fuel  gas  does  not  have 
sufficient  time  to  penetrate  in  depth  when  the  system  is  just 
started  up.  As  time  elapses,  however,  the  oxygen  diffuses  to 
the  regions  far  from  the  inlet.  Finally,  the  system  researches 


at  a  steady  state  after  experiences  a  time  period  of  transient 
variation. 

In  this  work,  the  surface  overpotential  0  represents  the 
voltage  loss  due  to  chemical  reaction  within  the  catalyst  layer. 
Therefore,  a  higher  surface  overpotential  means  a  fuel  cell 
with  a  faster  surface  reaction  rate.  To  explore  effects  of  the 
surface  overpotential  on  the  reactant  gas  transport,  the  con¬ 
tours  of  the  oxygen  mass  fraction  in  the  gas  diffuser  layer  and 
catalyst  layer  of  £\  =  0.5,  14-  0.225  and  A  =  0.75  at  different 
instants  are  shown  in  Fig.  3.  Comparison  of  Figs.  2  and  3 
discloses  that  a  smaller  oxygen  concentration  is  noted  for  a 
system  with  a  high  surface  overpotential  77.  This  is  resulted 
from  the  fact  that  a  PEM  fuel  cell  with  a  higher  surface  over¬ 
potential  represents  a  fuel  cell  with  a  faster  surface  reaction 
rate.  This  explains  why  the  dynamic  response  is  faster  in  the 
case  with  a  higher  77. 

To  examine  the  effects  of  the  channel  width  fraction  A 
on  the  reactant  gas  transport,  the  temporal  variation  of  the 
oxygen  mass  fraction  in  the  gas  diffuser  layer  and  catalyst 
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Fig.  4.  Oxygen  mass  fraction  contours  in  diffusion  and  catalyst  layers  with 
X  =  0.75  and  r)  =  0.2  at  different  transients:  (a)  f  =  0.01;  (b)  f  =  0.05;  (c) 
t  =  2.184. 


layer  of  s i  =0.5,  rj  =  0.2  and  A.  =  0.75  are  shown  in  Fig.  4. 
Comparing  with  the  results  in  Fig.  2,  it  is  shown  that,  in  the 
case  with  a  larger  channel  width  fraction,  the  distribution  of 
oxygen  is  more  uniform.  This  means  that  the  concentration 
gradient  of  the  oxygen  is  more  noticeable  in  the  case  of  a 
smaller  channel  width  fraction,  i.e.  A  =  0.5.  Additionally,  the 
time  needed  for  the  system  to  reach  the  steady  state  is  longer 
for  a  smaller  A.  This  can  be  made  plausible  by  noting  the  fact 
that  for  a  narrow  channel,  the  effective  inlet  region  is  smaller. 
In  this  situation,  it  takes  a  longer  penetration  time  to  diffuse 
oxygen  to  the  steady  state  when  the  system  is  started  up. 

Fig.  5  presents  the  local  oxygen  fraction  distributions 
along  the  interface  between  gas  diffuser  layer  and  catalyst 
layer,  i.e.  boundary  B4  in  Fig.  1,  under  various  conditions.  It 
is  noteworthy  that  in  the  early  state,  e.g.  t  =  0.01,  the  oxygen 
fraction  is  almost  zero.  This  means  that  the  oxygen  does  not 
penetrate  to  the  interface  at  this  short  duration.  Later,  the  oxy¬ 
gen  fraction  increases  with  time  till  steady  state  is  reached. 
With  the  periodic  conditions  at  left  and  right  sides,  i.e.  bound¬ 


0.03  - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - r 


Fig.  5.  Transient  variations  of  oxygen  mass  fraction  along  interface  between 
the  gas  diffuser  layer  and  catalyst  layer:  (a)  £1  =  0.5,  X  =  0.5  and  r /  =  0.2;  (b) 
£1  =0.6,  X  =  0.5  and  r\-  0.2;  (c)  £1  =0.5,  X  =  0.75  and  rj  =  0.2;  (d)  £1  =0.5, 
X  =  0.5  and  rj  =  0.225. 


aries  3  and  7,  the  distributions  of  the  mass  fraction  at  the  left 
sides  (i.e.,  v  <  0.5)  are  higher  than  that  at  right  sides.  A  careful 
examination  of  Fig.  5(a)  indicates  that  as  the  time  is  about 
0.194,  the  mass  fraction  is  very  close  to  that  at  steady  state 
(t  =  8.35).  It  is  found  from  a  separate  numerical  run  that  at 
t-  1.0,  the  mass  fraction  almost  collapses  with  that  at  steady 
state.  This  implies  that,  after  the  PEM  fuel  cell  system  is 
started  up,  it  arrives  at  the  steady  very  quickly.  Comparison 
of  the  corresponding  curves  of  £  =  0.01  s  in  Fig.  5(a)  and  (b) 
discloses  that  the  dynamic  response  to  the  power  surge  of  the 
PEM  fuel  cell  is  faster  in  the  initial  stage  for  a  system  with 
a  higher  GDL  porosity.  This  is  owing  to  the  fact  that  a  high 
GDL  porosity  enhances  the  diffusion  transport  of  the  reactant 
gas  through  the  porous  layers  and,  in  turn,  increases  the  local 
values  of  the  oxygen  mass  fraction.  However,  the  time  needed 
for  the  system  to  reach  the  steady  state  is  longer  for  a  PEM 
fuel  cell  with  a  higher  GDL  porosity.  As  explained  in  Fig.  3, 
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the  distributions  of  the  mass  fraction  along  the  interface  be¬ 
tween  the  GDL  and  catalyst  layer  are  relatively  uniform  for  a 
larger  channel  width  fraction,  see  Fig.  5(a)  and  (c).  Besides, 
a  shorter  dynamic  response  time  is  experienced  for  a  PEM 
fuel  cell  with  a  larger  channel  width  fraction.  Fig.  5(d)  shows 
the  transient  variations  of  the  mass  fraction  with  s i  =0.5, 
0  =  0.225  and  A  =  0.5.  Relative  to  Fig.  5(a),  it  is  observed  that 
an  increase  in  the  surface  overpotential  rj  results  in  a  decrease 
in  the  mass  fraction.  This  phenomenon  is  attributed  to  the  fact 
that  a  PEM  fuel  cell  with  a  higher  surface  overpotential  repre¬ 
sents  a  fuel  cell  with  a  faster  surface  reaction  rate.  Therefore, 
the  mass  fraction  along  the  interface  of  the  catalyst  layer  is 
lower.  In  addition,  the  response  time  is  shorter  for  a  PEM  fuel 
cell  with  a  higher  ij- 

In  Fig.  6,  the  unsteady  phase  potential  distributions  along 
the  interface  between  the  catalyst  layer  and  the  membrane  at 
the  conditions  corresponding  to  that  in  Fig.  5  are  presented. 
An  overall  inspection  on  Fig.  6  discloses  that,  in  the  early 
transient,  the  value  of  the  phase  potential  is  almost  zero,  i.e., 
the  value  of  initial  condition.  But  as  the  time  elapses,  the 
phase  potential  decreases  until  the  system  becomes  steady.  As 


Fig.  6.  Transient  variations  of  phase  potential  along  interface  between  the 
gas  diffuser  layer  and  catalyst  layer:  (a)  s\  =0.5,  Z  =  0.5  and  r\-  0.2;  (b) 
s i  =0.6,  A.  =  0.5  and  rj  =  0.2;  (c)  s\  =0.5,  X  =  0.75  and  *7  =  0.2;  (d)  s\  =0.5, 
X  =  0.5  and  rj  =  0.225. 


Fig.  7.  Transient  variations  of  local  current  density  along  interface  between 
the  gas  diffuser  layer  and  catalyst  layer:  (a)  s\  =  0.5,  X  =  0.5  and  r\  -  0.2;  (b) 
s i  =  0.6,  X  -  0.5  and  r\  -  0.2;  (c)  £\  -  0.5,  X  =  0.75  and  *7  =  0.2;  (d)  £1  =  0.5, 
X  =  0.5  and  *7  =  0.225. 

the  interpretation  for  the  results  in  Fig.  5,  a  shorter  dynamic 
response  time  is  noted  for  a  PEM  fuel  cell  with  a  higher  GDL 
porosity,  channel  width  fraction  and  surface  overpotential.  In 
fact,  with  the  results  of  the  phase  potential  distributions,  the 
average  phase  potential  difference  and  the  voltage  can  then 
be  obtained. 

The  distributions  of  the  local  current  density  along  the 
interface  at  various  instants  with  the  same  conditions  in 
Figs.  5  and  6  are  shown  in  Fig.  7.  Slow  response  of  the  sys¬ 
tem  in  the  early  initial  transient  is  again  demonstrated  by  the 
very  low  current  density.  With  the  time  going,  the  oxygen 
mass  fraction  penetrates  into  the  catalyst  layer  from  the  in¬ 
let.  Thereafter,  the  oxygen  reacts  and  produces  the  significant 
amount  of  current.  Additionally,  a  higher  current  density  is 
found  for  a  PEM  fuel  cell  with  a  higher  GDL  porosity,  channel 
width  fraction  and  surface  overpotential.  Since  the  increase  in 
GDL  porosity  enhances  the  diffusion  transport  of  the  reactant 
gas  through  the  porous  layers  and,  in  turn,  increases  the  lo¬ 
cal  values  of  the  current  density  with  a  specified  overpotential 
(by  comparing  Fig.  7(a)  and  (b)).  However,  this  high  transport 
rate  of  the  reactant  gas  in  the  across-the-cell  direction  may 
result  in  a  relatively  non-uniformity  in  x-distribution  of  the 
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Fig.  8.  Transient  variations  of  current  density:  (a)  A  =  0.5  and  ?7  =  0.2;  (b) 
£ i  =  0.5  and  r\  -  0.2;  (c)  £\  =  0.5  and  1  =  0.5. 

current  density  as  a  penalty.  A  careful  comparison  of  Fig.  7(a) 
and  (d)  indicates  that  the  distribution  of  the  current  density 
in  x-direction  is  relatively  uniform  at  low  overpotentials.  At 
a  high  overpotential,  the  surface  reaction  rate  becomes  fast 
and  the  relatively  weak  diffusion  leads  to  a  non-uniformity 
in  the  current  density  distribution. 

In  the  study  of  the  transient  behaviors  of  reactant  gas  trans¬ 
port  in  PEM  fuel  cell,  the  dynamic  response  time  for  the 
PEM  fuel  cell  subjected  to  a  power  surge  is  important  un¬ 
der  various  conditions.  Therefore,  the  temporal  variations  of 
the  current  density  and  voltage  for  different  conditions  are 
presented  in  Figs.  8  and  9,  respectively.  It  is  disclosed  that 
when  the  PEM  fuel  cell  is  subjected  to  a  power  surge  from 
the  stationary  condition,  the  PEM  fuel  cell  responds  quickly 
to  the  change.  Although,  the  full  time  duration  to  the  steady 
state  is  much  longer  than  1  s,  but  the  PEM  fuel  cell  almost 
changes  to  the  99%  value  at  steady  within  1  s.  In  addition, 
a  drastic  response  occurs  at  the  instant  of  about  0.1  s.  It  is 
found  from  the  separate  numerical  data  that  at  the  instant 
less  than  0.001  s,  the  PEM  fuel  cell  does  not  respond  to  the 
change  in  the  operating  condition.  This  means  that  the  PEM 
fuel  cell  is  still  stationary  with  almost  null  current  density  and 


Fig.  9.  Transient  variations  of  voltage:  (a)  A  =  0.5  and  rj  =  0.2;  (b)  s\  =0.5 
and  7)  =  0.2;  (c)  s\  =  0.5  and  1  =  0.5. 


voltage.  In  Figs.  8(a)  and  9(a),  the  effects  of  the  GDL  poros¬ 
ity  on  the  current  density  and  dynamic  response  time  can  be 
observed. 

In  the  study  of  the  transient  response  of  a  PEM  fuel  cell 
system,  the  time  constant  is  one  of  significant  parameters 
for  evaluation  of  the  system  performance,  especially  for  the 


Fig.  10.  Effects  of  GDL  porosity  s \  on  the  time  constant. 
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system  subjected  to  a  start-up  operation.  In  this  work,  the 
time  constant  is  defined  to  be  the  time  needed  for  the  current 
density  reaching  90%  of  the  value  at  the  corresponding  steady 
state.  Fig.  10  presents  the  effects  of  GDL  porosity,  channel 
width  fraction  and  surface  overpotential  on  the  time  constant. 
It  is  indicated  that,  in  the  parameter  range  considered,  the 
time  constant  is  less  than  0.4  s,  and  increases  in  GDL  porosity, 
channel  width  fraction  and/or  surface  overpotential  all  results 
in  a  reduction  in  time  constant. 


5.  Concluding  remarks 

The  dynamic  behavior  of  gas  transport  and  cell  perfor¬ 
mance  of  PEM  fuel  cell  has  been  examined  by  the  diffusion 
transport  model  developed  in  the  present  work.  With  the  as¬ 
sumption  of  the  two-dimensional  mass  transport  in  the  cath¬ 
ode  side  of  PEMFC,  the  effects  of  the  channel  width  fraction, 
A  =  fc/f  b,  the  porosity  of  the  gas  diffuser  layer,  £\ ,  and  the 
surface  overpotential  of  the  catalyst  layer,  77,  on  the  tran¬ 
sient  characteristics  of  the  resultant  current  density  and  mass 
transport  are  focused.  Based  on  the  results  and  analysis,  the 
following  conclusions  can  be  drawn. 

(1)  Surface  overpotential  r\  represents  the  voltage  loss  due  to 
chemical  reaction  within  the  catalyst  layer.  A  higher  sur¬ 
face  overpotential  means  a  fuel  cell  with  a  faster  surface 
reaction  rate.  Therefore,  a  faster  response  is  experienced 
for  the  PEM  fuel  cell  when  it  is  started  up  from  the  sta¬ 
tionary  condition. 

(2)  For  a  larger  channel  width  fraction,  a  faster  response  as 
well  as  better  cell  performance  can  be  attained. 

(3)  Morphology  of  the  porous  layers  such  as  the  gas  diffusion 
layer  in  the  present  half-cell  model  is  characterized  by  the 
porosity.  The  present  results  with  porosity  in  the  range 
of  0. 3-0.7  show  that  the  full  response  time  is  longer  for  a 
PEM  fuel  cell  with  a  higher  GDL  porosity.  While  the  time 
constant  is  shorter  for  a  cell  with  higher  GDL  porosity. 

(4)  It  is  disclosed  that  an  increase  in  A,  or  77  may  lead  to 
a  faster  dynamic  response  of  a  PEM  fuel  cell  subjected 
a  started-up  operation.  Although  the  dynamic  response 
time  may  be  as  long  as  10  s  due  to  the  mass  transport  lag, 
the  fuel  cell  system  needs  only  about  0.4  s  to  reach  the 
90%  of  the  current  density  at  the  corresponding  steady 
state. 


Acknowledgement 

The  study  was  supported  by  the  National  Science  Council, 
the  Republic  of  China  through  the  grant  numbers  NSC  92- 
2212-E-002-096,  NSC  91-2218-E-035-002,  NSC-92-2212- 
E-035-027  and  NSC  92-2623-7-002-006-ET. 


References 

[1]  D.N.  Bernadi,  M.W.  Verbrugge,  Mathematical  model  of  a  gas  dif¬ 
fusion  electrode  bonded  to  a  polymer  electrolyte,  AIChE  J.  37  (8) 
(1991)  1151-1163. 

[2]  G.  Xie,  T.  Okada,  Water  transport  behavior  in  Nation  117  mem¬ 
branes,  J.  Electrochem.  Soc.  142  (9)  (1995)  3057-3062. 

[3]  G.  Xie,  T.  Okada,  Pumping  effects  in  water  movement  accompanying 
cation  transport  across  Nation  117  membranes,  Electrochim.  Acta  41 
(9)  (1996)  1569-1571. 

[4]  T.  Okada,  G.  Xie,  M.  Meeg,  Simulation  for  water  management  in 
membranes  for  polymer  electrolyte  fuel  cells,  Electrochim.  Acta  43 
(14-15)  (1998)  2141-2155. 

[5]  T.  Okada,  G.  Xie,  O.  Gorseth,  S.  Kjelstrup,  N.  Nakamura,  T. 
Arimura,  Ion  and  water  transport  characteristics  of  Nation  mem¬ 
branes  as  electrolytes,  Electrochim.  Acta  43  (1998)  3741-3747. 

[6]  T.  Okada,  Theory  for  water  management  in  membranes  for  polymer 
electrolyte  fuel  cells.  Part  1.  The  effect  of  impurity  ions  at  the  an¬ 
ode  side  on  the  membrane  performances,  J.  Electroanal.  Chem.  465 
(1999)  1-17. 

[7]  T.  Okada,  Theory  for  water  management  in  membranes  for  polymer 
electrolyte  fuel  cells.  Part  2.  The  effect  of  impurity  ions  at  the  cath¬ 
ode  side  on  the  membrane  performances,  J.  Electroanal.  Chem.  465 
(1999)  18-29. 

[8]  V.  Gurau,  F.  Barbir,  H.  Liu,  An  analytical  solution  of  a  half-cell 
model  for  PEM  fuel  cells,  J.  Electrochem.  Soc.  147  (12)  (2000) 
4485-4493. 

[9]  S.  Dutta,  S.  Shimpalee,  J.W.V.  Zee,  Numerical  prediction  of  mass- 
exchange  between  cathode  and  anode  channels  in  a  PEM  fuel  cell, 
Int.  J.  Heat  Mass  Transfer  44  (2001)  2029-2042. 

[10]  J.S.  Yi,  T.V.  Nguyen,  Multi-component  transport  in  porous  electrodes 
of  proton  exchange  membrane  fuel  cells  using  the  interdigitated  gas 
distributors,  J.  Electrochem.  Soc.  146  (1)  (1999)  38-45. 

[11]  W.  He,  J.S.  Yi,  T.V.  Nguyen,  Two-phase  flow  model  of  the  cathode 
of  PEM  fuel  cells  using  interdigitated  flow  fields,  AIChE  J.  46  (10) 
(2000)  2053-2063. 

[12]  R.J.  Kee,  P.  Korada,  K.  Walters,  M.  Pavol,  A  Generalize  model  of  the 
flow  distribution  in  channel  networks  of  planar  fuel  cells,  J.  Power 
Sources  109  (2002)  148-159. 

[13]  E.  Hontabon,  M.J.  Escudero,  C.  Bautista,  PL.  Garcia- Ybarra,  L. 
Daza,  Optimisation  of  flow-field  in  polymer  electrolyte  membrane 
fuel  cells  using  computational  fluid  dynamics  techniques,  J.  Power 
Sources  86  (2001)  363-368. 

[14]  S.  Um,  C.Y.  Wang,  K.S.  Chen,  Computational  fluid  dynamics  mod¬ 
eling  of  proton  exchange  membrane  fuel  cells,  J.  Eelectrochem.  Soc. 
147  (12)  (2000)  4485-4493. 

[15]  W.M.  Yan,  C.Y.  Soong,  F.  Chen,  H.S.  Chu,  Effects  of  flow  distributor 
geometry  and  diffusion  layer  porosity  on  reactant  gas  transport  and 
performance  o  f  proton  exchange  membrane  fuel  cells,  J.  Power 
Sources  125  (2004)  27-39. 

[16]  M.  Francesco,  E.  Arato,  Start-up  analysis  for  automotive  PEM  fuel 
cell  systems,  J.  Power  Sources  108  (2002)  41-52. 

[17]  S.  Yerramalla,  A.  Davari,  A.  Feliachi,  T.  Biswas,  Modeling  and  sim¬ 
ulation  of  the  dynamic  behavior  of  a  polymer  electrolyte  fuel  cells, 
J.  Power  Sources  124  (2003)  104-113. 

[18]  S.  Um,  C.Y.  Wang,  Three-dimensional  analysis  of  transport  and 
electrochemical  reactions  in  polymer  electrolyte  fuel  cells,  J.  Power 
Sources  125  (2004)  40-50. 

[19]  F.  Chen,  Y.G.  Su,  C.Y.  Soong,  W.M.  Yan,  H.S.  Chu,  Transient  be¬ 
havior  of  water  transport  in  the  membrane  of  PEM  fuel  cell,  J. 
Electroanal.  Chem.  566  (2004)  85-93. 

[20]  S.V.  Patankar,  Numerical  Heat  Transfer  and  Fluid  Flow, 
Hemisphere/McGraw-Hill,  New  York,  1980,  Chapter  6. 


